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Experimental momentum distributions for the transitions 
to the ground state and first excited state of 6 He have been 
measured via the reaction 7 Li(e, e'p) 6 He, in the missing mo- 
mentum range from -70 to 260 MeV/c. They are compared 
to theoretical distributions calculated with mean-field wave 
functions and with variational Monte Carlo (VMC) wave 
functions which include strong state-dependent correlations 
in both 7 Li and 6 He. These VMC calculations provide a 
parameter- free prediction of the momentum distribution that 
reproduces the measured data, including its normalization. 
The deduced summed spectroscopic factor for the two transi- 
tions is 0.58±0.05, in perfect agreement with the VMC value 
of 0.60. This is the first successful comparison of experiment 
and ab initio theory for spectroscopic factors in p-shell nuclei. 

PACS numbers: 21.10.Jx, 21.60.Ka, 25.30. Dh, 27.20. +n 

In this letter we present new experimental data on the 
reaction 7 Li(e, e'p) 6 He and compare the deduced momen- 
tum distributions with recent ab initio predictions from 
variational Monte Carlo (VMC) calculations. Some years 
ago, two of us (Lapikas and Wesseling) participated in 
electron scattering experiments to determine shell oc- 



cupancies in the nuclei JU Si, 3 P, and J2 S [ljg]. In the 
latter Li2S target was used, with resolution suffi- 

cient to separate the discrete transitions in the reaction 
32 S(e, e'p) 31 P from those in the reaction 7 Li(e, e'p) 6 Hc. 
However, the results for 7 Li were not published at the 
time. Independently, one of us (Wiringa) recently calcu- 
lated the overlap wave function ( 6 He| a(r) | 7 Li) as part of 
a general program of quantum Monte Carlo studies of the 
light p-shell nuclei f|H]- These calculations use realistic 
two- and three-nucleon interactions fit to NN scattering 
data and few-body nuclear bound states, and produce 
fully correlated A-body wave functions. We have now 
found that using the VMC overlap as input to a Coulomb 
Distorted Wave Impulse Approximation (CDWIA) code 
results in excellent predictions of the observed momen- 
tum distributions and transition rms radii including the 
absolute normalization of the cross sections. This is the 
first successful comparison of experiment and ab initio 
theory for spectroscopic factors in p-shell nuclei. 

The effect of including short-range and tensor corre- 
lations in the calculation of nuclear structure has been 
studied previously in detail for few-body systems. In par- 
ticular, momentum distributions for the nuclei 2 H |^|, 
3 He 0, and 4 He H[|], measured via the reaction [e, e'p), 
have been compared to calculations (Faddeev [fill and 



VMC 0JT2|) that include state-dependent correlations 
derived from bare nucleon-nucleon interactions. Experi- 
ments on complex (A>4) nuclei have been performed E3] 
abundantly, but theoretical calculations (Green's func- 
tion [H and cluster VMC JT5|]) are more difficult, and 
have been limited to a few closed-shell nuclei; as discussed 
below, these usually overpredict the normalization of the 
cross sections. Typically, these data (consisting mainly 
of knockout of valence protons) are analyzed by compar- 
ing to calculations based on mean-field theory (MFT) 
that do not include (short-range and other) correlations, 
and by identifying the required renormalization as the 
spectroscopic factor. The resulting factors [jl3|,0 (about 
60-65% for nuclei ranging from A=6 to 209), are usually 
interpreted as evidence for the presence of important cor- 
relations. 

A quenching of this kind was predicted for infinite nu- 
clear matter [[17 18 but the extension of this result to 
finite nuclei is not straightforward due to the coupling 
to surface vibrations that affects the strength for transi- 
tions near the Fermi edge. For the nucleus 16 the effect 
of both short and long-range correlations was calculated 
with a Green's function method jl4j resulting in a reduc- 
tion of the strength to 0.76 of the MFT value. However, 
the inclusion of center-of-mass effects will probably in- 
crease this value to ~ 0.81 |lq], which is still consider- 
ably larger than the observed |1S|] lp strength (~0.6) at 
small excitation energies. It may be that 16 is an ex- 
ceptionally difficult case; more success has been gained 
with the larger nuclei 48 Ca and 90 Zr pOfl , although these 
calculations must use a G-matrix representation of the 
NN force, a step away from the bare interaction. In the 
present case, the VMC method uses the bare interactions 
to produce rather sophisticated six- and seven-body wave 
functions, including strong state-dependent correlations, 
which show the clustering expected in light p-shell nuclei. 

The theoretical description of the reaction (e, e'p) has 
been given in detail elsewhere |2l}| . In Plane- Wave Im- 
pulse Approximation (PWIA), the expression for the 
cross section reads: 
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(1) 



where the spectral function S(E m , p m ) denotes the prob- 
ability to find a proton with separation energy and mo- 
mentum (E m ,p m ) in the nucleus. The quantity Ko~ ep 
is the product of a phase space factor and the elemen- 
tary off-shell electron-proton cross section [p2l that de- 
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scribes the coupling of the virtual photon to the pro- 
ton. In this paper we consider the transitions to the 
0+ ground state (E m =9.97 MeV) and 2+ first excited 
state (£ ro =11.77 MeV) in 6 He. Hence in Eq. (1) E m 
has two discrete values and we can obtain the momen- 
tum distribution p(p m ) for each transition by integrating 
S(E m ,p m ) over the appropriate peak in E m . In PWIA 
the momentum distributions arc related to the overlap 
wave function ( 6 He| a(r) | 7 Li) via 

p(p m ) = | / e lp '"- r ( 6 He | o(r) | 7 Li)dr | 2 . (2) 

For the overlap wave functions we take either the MFT 
or the VMC results, which are discussed below. In order 
to account for Coulomb distortion of the electron wave 
functions and for final-state interaction (FSI) between 
the outgoing proton and the residual 6 He nucleus we use 
a CDWIA procedure p3| . Here the FSI are treated via 
an optical-model potential, the parameters of which were 
taken from a 100 MeV proton scattering experiment on 
6 Li p4j . For the extrapolation of these parameters to 
90 MeV protons and 6 He we used Schwandt's Q global 
parameterization for a large number of nuclei and ener- 
gies. The uncertainty due the treatment of the FSI was 
estimated by also extrapolating the optical-model param- 
eters from proton scattering data at lower and higher 
p6| , p7| energies. This yielded model uncertainties on the 
spectroscopic factors of 6% and on the deduced rms radii 
of the overlap wave functions of 2%. 
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FIG. 1. Momentum distributions in PWIA for the 
ground-state overlap wave function ( 6 He| o(r) | 7 Li) as calcu- 
lated in MFT (solid curve) and in VMC (circles). The error 
bars on the VMC data denote the uncertainty due to Monte 
Carlo sampling. The normalization of the MFT wave function 
has been chosen identical to that of the VMC one. 



In mean-field theory the overlap wave function 
( 6 He| a(r) Li) reduces to a single-particle wave function 
4> a { v ) since it is the product of two Slater determinants. 
It is taken as the solution of the Schrodinger equation 
with a Woods-Saxon potential that reproduces the ap- 
propriate binding energy. The radius of the potential is 
chosen such that the calculated momentum distribution 
fits the experimental data. Based on the VMC calcu- 
lations, which predict dominant and lpi/2 ampli- 
tudes for the two transitions, respectively, we choose a 
MFT IP3/2 wave function for the 3/2~ — > + ground 
state transition and a Xp\/2 MFT wave function for the 
3/2~ — » 2 + transition to the first excited state. 

Quantum Monte Carlo calculations of the ground and 
low-lying excited states for six- and seven-body nuclei 
have been made H using a realistic Hamiltonian con- 
taining the Argonne Vis two-nucleon |^9| and Urbana IX 
three-nucleon [£| potentials (AV18/UIX model). These 
calculations start with trial functions, J w , T), con- 
structed from products of two- and three-body correla- 
tion operators acting on a fully antisymmetrized set of 
one-body p-shell basis states that are LS coupled to the 
specified quantum numbers. Metropolis Monte Carlo in- 
tegration |3l| is used to evaluate (H>v I H | and 
diagonalize in the one-body p-shell basis, giving upper 
bounds to the energies of these states. The trial functions 
are then used as input to the Green's function Monte 
Carlo (GFMC) algorithm, which projects out excited 
state contamination in the trial function by means of the 
Euclidean propagation \1/(t) = exp[— (H— Eq)t]^v The 
GFMC results are believed to be within ~ 1% of the ex- 
act binding energy for the given Hamiltonian. 

For the AV18/UIX model, the GFMC energy for the 
7 Li ground state is -37.4(3) MeV, where the number in 
parentheses is the statistical error due to the Monte Carlo 
energy evaluation. The 6 He(0 + ) ground state is at - 
27.6(1) MeV and the 6 He(2+) excited state is at -25.8(2) 
MeV. While these ab initio energies are about 5% above 
experiment (which we attribute to inadequacies of the 
AV18/UIX model) the relative excitations of 9.8(3) and 
11.6(3) MeV are fairly close to experiment. The VMC en- 
ergies are not as good, but the one- and two-body VMC 
and GFMC density distributions are very similar, giving 
us some confidence in using the to study reactions. 
Recent calculations using the equivalent for 6 Li gave 
a very satisfactory description of both elastic and inelas- 
tic electron scattering form factors p2j. For the present 
work, the ( 6 He| <z(r) | 7 Li) overlaps were calculated using 
the techniques of Refs. 

In Figures |l| and || we compare the plane wave MFT 
and VMC wave functions in momentum space out to 
large momenta. In order to facilitate the comparison we 
have scaled the MFT overlap wave functions such that 
their normalization is identical to that of the VMC wave 
functions. For the ground state transition (see Fig. |l|) we 
observe that the MFT and VMC wave function are prac- 



2 



tically identical up to p m = 400 MeV/c, whereas above 
this momentum the VMC wave function is appreciably 
larger due to the inclusion of short-range and tensor cor- 
relations, which are absent in MFT. The transition to the 
first excited state contains both lp and If components, 
as shown in Fig. |2|. Here the deviation between the VMC 
and MFT wave function already starts at 250 MeV/c be- 
cause in MFT the wave function is purely lpx/2, whereas 
the VMC overlap contains four components (lpi/2, l.P3/2> 
1/5/21 1/7/2)- In addition to the effect of correlations 
these extra components cause an appreciable enhance- 
ment of the VMC wave function at high momentum rel- 
ative to the MFT wave function. 



7 Li(e,e'p) 6 He O VMC (1p+1f) 
3/2" ->2 + X VMC (1f) 




FIG. 2. Same as Fig. |l| but for the transition to the first 
excited state. The MFT wave function is pure lpi/2, the to- 
tal VMC wave function (circles) contains lp and If (crosses) 
components. 

The experiment was performed with the 1% duty 
factor electron beam from the NIKHEF medium- 
energy accelerator (MEA) and the high-resolution two- 
spectrometer setup in the EMIN end station [p3| . The 
data were taken concurrently with those for the reaction 
32 S(e,e'p) §§ for which purpose a self-supporting disc of 
Li 2 S was used as a target (thickness roughly 25 mg/cm 2 ). 
The target could withstand maximum average currents 
of 6 fiA when rotated continuously. The target thick- 
ness was monitored via frequent measurements of elastic 
scattering. The measurements were carried out in paral- 
lel kinematics for an outgoing proton energy of 90 MeV. 
As a result we needed two incident energies (329.7 and 
454.7 MeV) to cover the missing momentum range of -70 
to 260 MeV/c. Since the beam was tuned in dispersion 
matching mode |34| ] we could achieve an E m resolution 
of 180 keV (FWHM), sufficient to separate the discrete 



transitions from the two reactions. 

The data analysis was performed in a standard way de- 
scribed in detail elsewhere |3f| . From the measured cross 
sections we determined momentum distributions by in- 
tegrating over the appropriate missing-energy peak and 
by dividing out Ka epi for which we used the current- 
conserving expression cr^ 1 of De Forest ]22| . The result- 
ing experimental momentum distributions are shown in 
Fig. ||, where only the statistical errors are shown. The 
experimental systematic uncertainty on these data is 5%. 
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FIG. 3. Experimental momentum distributions for the 
transitions to the ground state (circles) and first excited state 
(crosses) in the reaction 7 Li(e, e'p) 6 He, compared to CDWIA 
calculations with MFT (solid) and VMC (dashed) wave func- 
tions. The dash-dot-dot curve represents the If contribution 
to the full VMC curve for the transition to the 2 + state. The 
error bars on the data are statistical only. For clarity data 
and curves for the ground state transition have been scaled 
by 10. 

In the only one earlier reported |36| study of the re- 
action 7 Li(e, e'p) 6 He the missing-energy resolution of 7 
MeV was insufficient to separate the two transitions pre- 
sented in this Letter. However, when corrected for the 

TABLE I. Spectroscopic factors (S) and rms-radii deduced 
in the present experiment for the transitions to the + and 2 + 
state in 6 He (first row). The listed errors include statistical, 
systematic and model uncertainties. The second (third) row 
presents the corresponding values for the VMC calculation 
with lp (lp+lf) wave function components. 



Model 


S 


S 


S 


rms [fm] 


rms [fm] 




0+ 


2+ 


0++2+ 


0+ 


2+ 


Exp. (lp) 


0.42(4) 


0.16(2) 


0.58(5) 


3.17(6) 


3.47(9) 


VMC (lp) 


0.41 


0.18 


0.59 


3.16 


3.14 


VMC (lp -+ 


•1/) 0.41 


0.19 


0.60 


3.16 


3.16 
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presence of some unresolved Is knockout strength and 
the difference in ejected proton energy, their momentum 
distribution, integrated over the region E m =6-15 MeV, 
agrees within error bars with that for the sum of the two 
transitions studied here. 

In order to compare the theoretical calculations with 
the data we carried out CDWIA calculations with the 
MFT and VMC wave functions as input. For the mean- 
field calculations we treated the normalization, i.e., the 
spectroscopic factor S, and the radius of the WS poten- 
tial (that fixes the rms radius < r 2 > 1 / 2 of the wave 
function) as free parameters to be determined from a 
least squares fit to the data. The resulting values are 
listed in table 1. The summed spectroscopic strength for 
lp knockout is 0.58±0.05, where we have included the 
experimental systematic uncertainty and the uncertainty 
due to the choice of the optical potential. The observed 
reduction of the single particle strength to 58% of the 
MFT value (which is unity for a single proton in the lp 
shell) is in good agreement with the reduction found for 
a large number of other complex nuclei [ji"3"f . 

Figure ^ also shows the calculated momentum distri- 
butions with the VMC wave functions, which are essen- 
tially parameter free. The agreement with the data is 
very good as shown in table 1 where the calculated spec- 
troscopic factors with these wave functions are given. 
The summed strength (0.60) for both transitions agrees 
within error bars perfectly with the value 0.58±0.05 de- 
duced from the MFT analysis. 

The VMC rms radius for the ground state transition 
agrees with the value deduced from the MFT analysis, 
showing that the calculated VMC ground state wave 
functions for 6 He and 7 Li have the correct shape. For 
the transition to the first excited state the rms radius of 
the VMC wave function is smaller than that found in the 
MFT analysis. This is caused by the different structure 
for both overlaps: the MFT wave function was assumed 
to be pure lpi/2, whereas the VMC wave function con- 
tains 1^3/2, I/5/2 an d I/7/2 components in addition. The 
contributions of the 1/ components, which depend sen- 
sitively on the details of the nucleon-nuclcon interaction 
employed, would show up in measurements at higher p m 
and could thereby serve as a further accurate test of the 
VMC wave functions. 

In summary we conclude that for the first time struc- 
ture calculations for a complex nucleus, based on a re- 
alistic nucleon-nucleon force, have been performed and 
compared to (new) experimental data for the reaction 
7 Li(e,e'p). The calculated spectroscopic strength (0.60) 
explains the reduction of the strength to 0.58±0.05 found 
in a MFT analysis of the data, while the calculated shape 
of the momentum distributions for lp transitions nearly 
coincides with the experimental data. Thus we have con- 
firmed the necessity of including full correlations in the 
nuclear wave functions. 
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